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Abstract 

Kinetic and magnetic enstrophy are shown to concentrate to smaller scales from the 
integral scale to a Kraichnan-type scale for the 3D magnetohydrodynamic equations. 
This is an improvement of the result from Bradshaw and Grujic (2013), using redesigned 
ensemble averages. 


1 Introduction 

Observational and numerical evidence suggests that in magnetohydrodynamic (MHD) tur¬ 
bulence the vorticity u = V x u and current j = V x b (here, u and b are the velocity and 
the magnetic held, respectively) concentrate on coherent quasi low-dimensional structures- 
predominantly quasi two-dimensional sheets-which become increasingly thin exhibiting mor¬ 
phological dynamics consistent with a process known as turbulent cascade. Understanding 
this process is important due to its intimate relationship with violent reconnection events 
in the realm of solar wind turbulence, and is highly relevant to the goals outlined in the 
Space Studies Board of the National Research Council (NRC) survey, 22A Decadal Strategy 
for Solar and Space Physics (Heliophysics), completed in 2012, and in particular, to the 
strategic goal to “discover and characterize fundamental processes that occur both within 
the heliosphere and throughout the universe” [17]. 

As a matter of fact-in the last hve years-there has been a flurry of activity in the 
heliophysics community directed at understanding the phenomenon of turbulent dissipation, 
especially within the range of kinetic scales. One of the most promising theories, supported 
by extensive computational simulations, is that the coherent structures (most notably current 
sheets) exhibit a process of turbulent cascade down to the kinetic scales, essentially, all the 
way to electron scales, where they trigger extremely strong and localized heating of the 
plasma, dissipating the energy (the so-called heating via current sheets) [14, 22, 23]. There 
is also a sense of optimism in part of the community that the relevance of this theory to 
the solar wind could be conhrmed via the data to be collected by the upcoming NASA 
Magnetospheric Multiscale (MMS) mission. 

The zero-step in this theory is an assumption on the geometry of the turbulent plasma 
at the interface between the continuum (described by the MHD system) and the kinetic (de¬ 
scribed by the Vlasov-Maxwell-Poisson system)-scale dynamics; shortly, the predominance 
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of the current sheet geometry originating at continuum scales is assumed, and then utilized 
as the input when descending into the kinetic-scale dynamics. This brings out the question 
of existence of turbulent cascades-and in particular, the cascades of kinetic and magnetic 
enstrophies-in the continuum/MHD regime to the forefront of scientihc interest. 

Kinetic and magnetic enstrophy was shown to concentrate towards smaller scales in [1] 
using ensemble averages in physical scales. Here the ensemble averages used are redesigned to 
allow weaker assumptions. The ensemble averages used to state the enstrophy concentration 
theorem are described hrst. The following section goes over the conditions under which 
enstrophy concentration is shown. Then we recall the bonnds obtained in [1], and hnally the 
enstrophy concentration theorem is proven. 

Related work on cascades and locality in hydrodynamic tnrbulence can be fonnd in 
[2-5,7-13,15,16,18-21]. 


2 Fluxes and ensemble averages 

The MHD equations, which model evolntion of the velocity and magnetic fields in an elec¬ 
trically conducting incompressible flnid, read 

dtu — Au + {u ■ V)u — {b ■ V)b + VP = 0, 

dtb-Ab+{u-V)b-{b-V)u = 0, (1) 

V ■ u = V ■ b = 0. 


(Here, P is the total pressnre, and the magnetic resistivity and the kinematic viscosity are 
normalized to 1.) 

Taking the cnrl of eqnation 1 gives eqnations for vorticity and cnrrent, 
dtco — Au = —{u ■ V)u {u ■ V)u + {h ■ V)j — {j ■ V)b, 

^ ( 2 ) 

dtj — Aj = —{u ■ V)j + {j ■ V)u + {b ■ V)u — {u ■ V)b + 2 ^ Vbi x Vui. 

i=i 


These eqnations can be nsed to stndy inward kinetic and magnetic enstrophy fluxes. 
Rather than traditional flnxes throngh the bonndary of a ball B = B{xq, 2R), 


/ -\u\^{u ■ n) da = 

'dB 2 

[ hj\‘^{u-n)da = ■ 

'dB ^ 


{u ■ V)a; ■ u dx, 


IB 


{u ■V)j ■ j dx, 


(3) 


'B 


we use a smooth cutoff function -ip supported on B{xo, 2R) and eqnal to 1 on B{xo, R)- 
with inward pointing gradient-and stndy fluxes through the spherical layer of thickness R, 


-\u\'^{u ■ Vip) dx = — J {u ■ V)a; • ipu dx, 
• V^/J) dx = - f {u- V)j ■ tpjdx. 


(4) 
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The reason is that this form of the flux is more amenable to mathematical analysis, while 
at the same time preserving the physics. Time-averaged quantities are studied in turbulence, 
so we take time averages weighted according to a smooth function r]{t) which is 0 on [0, T/3] 
and 1 on [2T/3,T], denoting (j){x,t) = 


-\u:\^{u -Vcf)) dx dt = — 


1 


\i\^{u ■ V0) dxdt = — 


{u ■ V)a; ■ (poj dx dt, 


{u ■ V)j • (f)j dx dt. 


(5) 


Multiplying equation 2 by (pu and (pj respectively, and integrating over space and time, 
will yield expressions which can be used to dynamically estimate the fluxes: 


-\u\‘^{u ■ Vp) dx dt = / -\u{x,T)\‘^^jJ{x) dx + 


I Va;p0 dx dt 


\u\^{ds(t) + Ap) dx dt — 


{b ■ V)j • {(pea) dx dt + 


(ca • V)m • {poo) dx dt 


{j ■ V)6 ■ {poo) dx dt 


( 6 ) 


\oo{x., T)\^p{x) dx + 


\Vu\^p dx dt + H^ + Np + L‘^ + Np 


./o 


■\j\‘^{u-Vp)dxdt= / -\j{x,T)\‘^p{x)dx + 


I dx dt 


■\jp{dsp + Ap) dx dt + 


{b ■ V)a; • {pj) dx dt — 

3 


-/ /(2Ev.,xV., 


= / ^\j{x,T)\‘^P{x)dx + 


J {oo ■ V)6 ■ {pj) dx dt 
{j ■ V)m ■ {pj) dx dt 

Jo J 

■ {pj) dx dt 

(*T p 

/ I Vjf 0 dx dt + W + Ni + U + Ni+X 


(7) 


Enstrophy concentration will be demonstrated locally, over a ball i?(0,2i?o)- Ultimately 
we want to show that the average of these fluxes over suitable collections of functions tp of 
a particular scale is positive, for a range of scales. 

Fix Uq > 1 and 3/4<p<l. A refined test function at scale R is any smooth function 
p supported in a ball of radius 2R satisfying 0 < 0 < 1, IV^I < ^0^, and |A-0| < 

Now £x a scale i?o refined test function pQ centered at 0. An ensemble at scale R 
with global multiplicity Ki and local multiplicity K 2 is a collection of scale R test functions 
{pijP^i satisfying the following properties: 


1 . pi <po <Y.Pi 

2. {Ro/R)^ <n< Ki{Ri^/R)^ 
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3. No point of B{2Rq, 0) is contained in more than K 2 of the supports of ipi. 

For a function /, denote by (-F)i? the ensemble average f 

Fo = ^ I /^o dx. 

Property 1 above is needed to compare (-F)_r to Fq. Due to Property 1, test functions 
near the boundary of the support of ipQ will have small integrals, effectively skewing the 
ensemble average towards zero. Larger Ki and K 2 allow ensembles that have higher weight 
on functions away from the boundary, making the skewing insignihcant. 

These ensemble averages can be viewed as a way to detect whether a function is signif¬ 
icantly negative above some spatial scale. If every ensemble average {F)ji is positive, no 
matter how one arranges and stacks the test functions, then the function is not signihcantly 
negative at scales larger than R. Increasing Ki and K 2 lowers the threshold for a function 
to be considered signihcantly negative. 

One way to explicitly construct ensembles is to apply Lemma 2 to 'ipo and varying the 
multiplicity of the resulting functions (Assume ipQ satishes the stronger Og-bounds to get an 
ensemble with Og-bounds). 

The following lemma states that ensemble averages (at any scale) of positive functions 
are comparable to the large scale mean. The proof immediately follows from the dehnitions. 

Lemma 1. /// > 0 then -^Fq < {F)r < K 2 FQ. For slightly modified ensemble averages, 
we have wf - ^ 2 ;^ / f'^o dx {S>0). □ 

Using a rehned partition of unity, one can turn larger scale ensembles into smaller scale 
ensembles. 

Lemma 2. Any seale R test function satisfying Cq bounds is a sum of 8 \R/R'~\^ scale R! 
test functions satisfying Cq bounds (where R> R! and C( depends only on Cq). 

Therefore for all {Ki, K 2 ,Co)-ensembles at scale R and every R' < R, there exists a 
{fi)AKi, 8 K 2 ,C'Q)-ensemble at scale R' such that {F)r = {F)r'. 

Proof. Let be a scale R test function satisfying Cq bounds . Now to construct the partition 
of unity, take a scale R' test function go (satisfying Oq bounds), centered at zero and equal 
to 1 on [—R', R']^. Dehne gp = go{x — 2R'p), where p G Ifi. Then 1 < J2p9p < 2 so we may 
dehne hp = gp/ Y.q 9 q- 

Some calculus shows that |V/ip| < and |Ahp| < so \'V{'fihp)\ < 

'^{jphpY and |A('^hp)| < ('^hp)^^~h Fewer than 8 \R/R'Y < 64(i?/i?')^ of the 

functions f)hp are nonzero, and for any x, 'fip{x) Y 0 most 8 functions. 

Since ip = '^pfi’hp, the hrst claim is proven. For the second claim, given an ensemble 
{ipi}i, the new ensemble will be {ipihp}i^p. □ 


3 Assumptions 

Now we come to the conditions on the current and vorticity over B{0,2Rq) x (0,T) under 
which we can show that there is enstrophy concentration. 
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3.1 Geometry/smoothness 

Denote by 9{u, v) the angle between two vectors u, v. It is required that for some M, Ci > 0, 

I sir).e{uj{x + y,t),uj{x,t))\ < Ci\y\^/‘^ (8) 

for every t G (0,T), every x G 5(0, 2i?o + with \'Vu{x,t)\ > M, and every y such 
that \y\ < 2((To// 9) + (o'o//9)^'^^, and 

\i{x + y,t)- j{x,t)\ < \i{x + y,t)\\y\^/‘^ (9) 

for every t G (0,T), every x G 5(0,2i?o + -Rq'^^) with \'Vb{x,t)\ > M, and every y such that 
\y\ < 2((To// 9) + ((Jo//9)^'^^. The quantity Uq/(3 is defined in section 3.2. 

Condition 8 depends only on the angle of the vorticity vector u. This ^-Holder coherence 
will deplete the vortex stretching term. This condition is needed because condition 9 is 
insufficient to control the vortex stretching term which has no explicit dependence on 
the magnetic held. 

Condition 9 is needed for the nonlinear terms that do not have a geometric kernel avail¬ 
able. It requires ^-Holder continuity, and more. It would be too restrictive if \j{x + y,t)\ 
were ever small, but x + y is always close to the region where | V6| is large, and roughly, V6 
and j are large in the same regions. 

In [1], condition 9 was used with j replaced by u with no assumption corresponding to 
8. The present formulation is preferred because current appears to be more regular than 
vorticity in numerical simulations. 

3.2 Kraichnan-type scale 

Let Co, Eq, and Pq denote the time-averaged total energy, total enstrophy, and total palin- 
strophy at the integral scale. Precisely, 



Dehne the modihed Kraichnan-type scale (Tq by 



Assumption 2 is that cq < /3Ro, where /? is a constant (0 < /9 < 1) identihed in the 
proof. If u and j have large gradients as expected in a turbulent how, this assumption will 
be satished. 
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3.3 Localization 

Any smooth solution to the equations will have ca G L^((0,T) x 5(0, 2i?o + It is 

required that the kinetic and magnetic enstrophy is not too highly concentrated around any 
point: 



+ IjP dx dt < 


1 


( 11 ) 


for any y G 5(0, 2Rq) where Ci is a constant and R = 2aojf3 + ((To// 9)^/^. 
In [1], it was required that 


/o JB{ 0 , 2 Ro+rI/^) 


U) 


P dx dt < —. 


This assumption restricts the extent of the range of scales over which we can show enstro¬ 
phy concentration. The magnetic enstrophy is in the present assumption because of the 
geometric/smoothness assumption on the current that is not found in [1], 


3.4 Modulation 


The assumption imposes a restriction on the time evolution of the integral-scale kinetic 
and magnetic enstrophies across (0,T) consistent with our choice of the temporal cutoff. 
Precisely, 


\uj{x, T)l^'^o(a^) dx > - sup 

2 t . 


\j{x,T)\‘^ijo{x) dx> - sup 

^ t , 


\uj{x^ t)\^il)Q{x) dx, 
|j(a;,t)|Vo(a:) dx, 


4 Bounds 


In [1], the terms of equations 6 and 7 are bounded by quantities which can be related to 
eo,5o, and Pq. 

Obtaining the desired bounds using Assumption 3.1 will use all of the same techniques as 
in [1] with the major difference being labelling, except for the vortex stretching term which 
uses geometric depletion as in [6]. 

Ultimately, we have 


5"^ + W + + N( + iV| + L‘^ + L^ + X 



j '4){x){\u{x,t)\^+ \j{x,t)\^)dx + 
-|- |jp) dx dt + 


P4 



04P-3 


Mp -I- l&P 
2 


dx dt. 
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where Kp,Ke, are constants and a is an interpolation parameter that we may choose. 
The constants Kp, Kp, K^. are fixed upon a choice of the parameters Ki, K 2 , p,Cq,Ci, and 
M. In the following section Kp, Kp, are instead the constants obtained corresponding to 
64i^i,8i^2, and Cq. 


5 Main result 


Kinetic and magnetic enstrophy is on average being transported to smaller scales down to the 
Kraichnan-type scale (Tq/13. The following notation will be used in the proof: For a density 
/ and ensemble the average over one element of the ensemble is Fi = J fipidx and 

the ensemble average is {F)p = Sr=i ■ 00 + {u ■ V)j ■ j dt, 

the enstrophy flux density, so that = J (pipi dx is the kinetic and magnetic enstrophy flux 
centered at Xi at scale R. Referring to equations 6 and 7, suppressing the subscript i, denote 
H = H‘^ + and likewise for N and L. Define similarly 


and 


e = ^ / + \b\^) dt dt, 

T 

P = j ^lJ,{x){\u{x,T)\^ + \J{x,T)\^)dx + j (j),{\Vu\^ + \Vj\^)dxdt, 


P = 1 sup f ^jJi{x){\u{x,t)\^+ \j{x,t)\‘^)dx + f f (l)i{\Vu\^+\Vj\‘^)dxdt. 

^ te(o,T) J Jo J 

Theorem. Under Assumptions 1-4, for any Ki and K 2 there exists such that for any 
{Ki, K 2 )-ensemble at scale R ranging from a^/(3 to Rq, we have —Po < {^)r < K*Po- 

Proof. We start by showing that for R = (Jo// 9, for (64iFi, 8 iF 2 )-ensembles satisfying Cq 
bounds, we have —Po < {^)r < K*Po- 

Referring to equations 6 and 7, ( 4 *)/? = {P)r + {H + N + L + X)r, where p is the positive 
density ||a;(T)p + \Wuj\^'qdt so that q^Po < {P)r < 8K2P0. The rest of the terms are 

relatively small upon averaging: |Tr + iV + L + X| < Kp(^ — h \ \uj\\L2(B{xi,2R+R'^/^)x{o,T))^P + 
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—E + -^e, so that 


\{H + N + L + A')h| < (A>(i + M\)P + ^E + ^e) 


R 


< 8 K 2 Kp[ —h ||a;|| ]Pq + 

.a 

1 


, a^8K2Ke 
Eo + --eo 


i ?2 ^ ' i?4 

= (8K2Kp{^ + ||a;||) + 8K2KeI3^ + a^8K2K^(3^')Pi 


< 


A-UKi 


-Pn 


by taking a = 4 ■ 64Ki ■ 8 K 2 EP, using that (4 is sufficiently small such that 8 K 2 KeI 4 ‘^ + 

■ 8 K 2 Ke/ 3 ^ < and using the assumption that ||w||L 2 (B(^,, 2 ij+ij 2 / 3 )x(o,T)) < ^ = 

1 


4-64iCi-8iC2i^F ■ 

Therefore J Pq < {^)r < { 8 K 2 + ^ )Po for any (64ifi, 8 K 2 , Co)-ensemble at 

4 ■ d4Ai 4 ■ d4Ai 

1 

scale R = (Jq/(3. For the range of scales from a^/(3 up to Rq^ by Lemma 2, - t-ttPo < 

4 ■ 64Ai 

3 

{^)r < { 8 K 2 + - -7— )Po for all scale R (iFi, 7 ^ 2 , Co)-ensembles. 

4 ■ 64ai 


□ 


The locality of magnetic and kinetic enstrophy flux stated in [1] also holds in our setting 
of more satisfying assumptions and modihed ensemble averages. The result is stated in terms 
of the time-averaged enstrophy flux, rather than the time-averaged enstrophy flux per unit 
mass used above. That is, 


for an ensemble at scale R. 

Corollary. Under Assumptions 1-4, for any Ki,K 2 there exists K^, such that for any r,R 
between Co/fl and Rq and any {Ki,K 2 ) ensembles, enstrophy flux is local: 


1 / r 

'i^Sr 


< 


{^)r 

{'^)r 


< k: 
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